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i.  introouctiou 

Zn  1933  Kapitza  and  Dirac  [1]  predicted  that  electrons  traveling  in  a 
standing  light  wave  would  be  reflected  from  the  planes  of  peak  intensity.  The 
probability  per  electron  for  reflection  was  proportional  to  the  square  of  the 
product  of  the  field  intensity  and  the  interaction  tine.  Because  the  light 
sources  available  at  that  tine  were  too  weak  to  generate  an  observable  effect, 
the  subject  was  neglected  until  the  discovery  of  the  laser.  Since  then, 
numerous  theoretical  and  experimental  studies  have  been  devoted'  to  this  subject. 

The  earliest  of  these  theoretical  papers  [2,3]  followed  closely  the 
original  treatment  of  Kapitza  and  Dirac  inasmuch  as  they  relied  on  a  first-order 
expansion  of  the  wavefunction  in  terms  of  the  standing  wave  field.  Like  the 
original  paper,  they  are  not  valid  for  presently  available  laser  intensities 
sinoe  they  would  predict  scattering  probabilities  in  excess  of  unity.  M.  V. 
Fedorov  14]  published  the  first  extensive  treatment  of  the  problem.  By 
neglecting  the  time  dependent  part  of  the  standing  wave,  he  was  able  to  rewrite 
the  Schrodinger  equation  in  the  form  of  a  Mathieu  equation.  Solutions  were  then 
found  for  the  cases  of  a  low  intensity  field  and  a  high  intensity  field. 
Unfortunately,  neither  of  these  cases  corresponded  to  intensities  used  in  the 
experiments. 

Gush  and  Gush  [5]  used  the  nonrelativistic  Green's  function  for  an  electron 
in  a  standing  wave  field  to  produce  an  exact  solution  to  the  problem  when  the 
time  dependent  part  of  the  field  is  neglected.  This  treatment  is  not  only  valid 
for  all  intensities  up  to  the  point  where  the  time  dependent  part  of  the  field 
becomes  important  but  also  for  electron  momenta  that  do  not  satisfy  the  Bragg 
conditions.  Furthermore,  the  probabilities  for  higher-order  reflections  are 
treated.  Unfortunately,  the  final  expressions  for  the  scattering  probabilities 
are  not  in  terms  of  known  functions  and  are  unwieldy. 
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The  question  of  what  role  if  any  the  tine  dependent  part  of  the 
standing-wave  field  plays  in  the  scattering  process  is  still  an  unresolved 
issue.  Gush  and  Gush  [5]  and  others  have  argued  that  the  tine  dependent  portion 
of  the  field  can  safely  be  neglected  for  the  intensities  used  in  the  reported 
experiments.  This  conclusion  has  been  disputed  by  Ehlotzky  et  al.  [6]  whose 
calculations  show  that  the  tine  dependence  of  the  field  has  a  significant 
influence  on  the  scattering  of  the  electron.  The  validity  of  their  calculation 
remains  in  doubt  in  view  of  the  fact  that  except  for  snail  interaction  times  and 
intensities  the  scattering  probabilities  for  the  time  averaged  case  differ  from 
those  predicted  by  the  more  exact  treatment  of  Gush  and  Gush  15] . 

In  this  paper,  the  role  of  the  time  dependent  portion  of  the  standing  wave 

is  investigated  using  several  different  approximation  techniques.  The  geometry 

and  governing  equations  are  introduced  in  Section  II.  Various  nondimensional 

parameters  are  also  defined  in  this  section  and  interpreted  in  terms  of  physical 

quantities.  In  Section  III  multiple  time-scale  perturbation  theory  is  used  to 

calculate  scattering  probabilities  for  the  case  when  the  standing  wave  field  is 

not  too  strong.  This  case  corresponds  to  the  one  previously  treated  in  the 

literature.  The  time  averaged  equation  is  derived  with  no  further  assumption 

and  the  limits  of  its  validity  are  discussed.  This  equation  is  then  solved  in 

terms  of  Mathieu  functions  and  several  new  features  of  the  scattering 

probabilities  are  discussed.  In  Section  IV,  we  treat  the  case  that  is 

characterized  by  strong  coupling  and  high  electron  energies.  This  case  is 

discussed  using  a  semiclassical  approximation.  The  wavefunction  is  again  found 

to  be  quasiperiodic  within  the  limitations  of  perturbation  theory.  Under 

appropriate  conditions  the  corresponding  classical  problem  leads  to  focusing. 

In  quantum  mechanical  terms  the  focusing  is  exhibited  as  sharp  localized  maxima 

is  the  scattering  probabilities.  In  Section  V  the  conclusions  are  drawn.  , ,  - , 
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II.  MODEL 


In  this  article  we  will  consider  the  problem  of  a  nonrelativistic, 
quantum-mechanical  electron  interacting  with  a  classical,  electromagnetic, 
standing  wave  as  illustrated  in  Figure  1.  The  electron  has  initial  energy  E  and 
momentum  p  which  is  chosen  so  as  to  lie  in  the  x'-z'  plane.  A  standing  wave  of 
frequency  u  lies  along  the  xVaxis  and  is  assumed  to  be  plane  polarized  along  the 
y'-axis.  The  interaction  is  turned  on  at  t  ■  0  and  continues  for  a  time  T  «  L/v 
where  v  «  |p|/m  is  the  electron  velocity,  m  is  the  electron  mass,  and  L  is  the 
interaction  length. 

The  standing  wave  is  given  by  the  vector  potential, 

X  *  A(cos(kx*-a)t)  +  cos(kx’  +  ut))y  *  2A  coskx'  cosut  y.  (2.1) 

Here  A  is  the  intensity  of  one  of  the  two  counter-propagating  fields  that 
combine  to  make  the  standing  wave  and  k  *  u/c» 

Because  the  electron  is  nonrelativistic  and  spin  effects  are  not  important 
for  most  cases  of  interest,  its  behavior  is  described  by  the  solution  of  the 
Schrodinger  equation  with  the  external  field  (2.1).  This  equation  can  be 
written  in  the  following  conventional  nondimensional  form: 

32  3Ui 

I -  -  2q  cos2x (1  +  cos2Ve)l  ♦  *  -  i  TT  (2.2a) 

ax2 


where 


For  all  cases  in  which  an  experiment  is  feasible  «  me2  so  that 


e  «  1  (2.3) 

The  initial  condition  for  the  wave  function  is  specified  by  its  free-spaoe  value 
at  the  time  when  the  interaction  is  turned  on 

#(x,t  «  0)  ■  exp(iBx)  (2.4) 

where 


6  -  Px/hk  (2.5) 

A  cursory  inspection  of  Eq.  (2.2a)  shows  that  two  types  of  interaction  of 
the  electron  with  the  field  can  occur.  In  one  case  the  particle  only  exchanges 
momentum  with  the  field.  This  is  the  elastic  scattering  mode  that  has  been 
extensively  investigated  11-5].  The  term  in  the  potential  that  only  has  a 
spatial  dependence  determines  this  interaction.  The  remaining  term  is  a  product 


of  a  spatial  and  temporal  piece.  Here  the  electron  experiences  a  change  in  both 
momentum  and  energy  each  time  it  interacts  with  the  field.  This  is  the 
inelastic  portion  of  the  scattering. 

Because  the  potential  in  Eq.  (2.2a)  only  permits  the  electron  to  change  its 
nondimensional  momentum  by  multiples  of  2,  the  final  state  of  the  electron  must 
be  characterized  by  a  momentum  pn  which  satisfies 


pn  ■  hkBn  ■  hk(B  +  2n),  n  ■  0,  ±1,  ±2,  ...  (2.6) 

The  probability  amplitude  that  after  an  interaction  time  rf  the  electron  is  in  a 
state  with  momentum  pn  is  given  by 


Pn<6'V 


2 

exp  (-i  6nz)  th  (z ,  x  f )  dz 


(2.7) 


These  probabilities  must  satisfy  the  condition 


l  P„  -  1  (2.8) 

n*-«p 


The  strength  of  the  interaction  is  determined  by  the  parameter  q  introduced 
in  Eq.  (2. 2d).  In  order  to  provide  a  physical  interpretation  for  this  quantity 
it  is  convenient  to  rewrite  it  in  terms  of  the  individual  wave  intensities 
I  »  ck2A2/8ir  so  that 


q 


2irr0l2I 

Kwc 


2mc2 

ftu 


(2.9) 


where  r0  is  the  classical  electron  radius  and  X  *  1A«  When  I  is  rewritten  in 
terms  of  the  photon  number  density  p  (I  *  hu>cp)  then  the  quantity  in  the  square 
brackets  reduoes  to  2nrnX2p.  The  coupling  constant  q  is  therefore  proportional 


to  the  nuriber  of  photons  in  a  cylinder  of  length  2ro  and  radius  1.  For 
available  laser  sources  q  can  reach  values  of  ^10 8  in  unfocussed  beams  although 
q  *10  is  tribre  characteristic  of  the  published  experiments  [7-12]. 

Because  we  have  been  unable  to  find  an  exact  solution  for  Lg.  (2.2a) ,  we 
have  had  to  resort  to  various  perturbation  schemes  to  find  approximate 
expressions  for  the  wave  function  in  different  regimes.  These  regimes  are 
characterized  by  the  relative  size  of  the  adjustable  parameters  q,  e  and  6.  In 
the  remainder  of  this  paper  we  will  be  concerned  with  the  range  of  values  of  the 
parameters  accessible  to  the  experimentalist. 

111.  WEAK  COUPLING 

The  weak  coupling  limit  is  characterized  by  the  two  conditions 

1/e  »  q 
and 

1/e  »  |el. 

This  case  covers  the  conditions  that  characterize  previous  experiments  [7-12] , 
i.e.  |6|  <10;  q  <  10  and  e_1,*106.  Under  these  conditions  multiple-time-scale 
perturbation  theory  [13]  can  be  used  to  find  an  approximate  expression  for  the 
wavef unction. 

In  order  to  apply  this  perturbation  analysis,  the  wavefunction  ip  is  assumed 
to  have  the  following  asymptotic  expansion: 

N 

Mz,t)  -  l  enip (x,t  )  +  0(eN+1)  ($.1) 

n«0 

where  each  ip  is  a  function  of  the  time  scales 


The  reason  for  considering  a  sequence  of  time  scales  is  that  each  mode  of  the 
wavefunction  has  a  characteristic  time  scale  £  ■  o>  (e ) t  that  depends  on  e.  The 
most  convenient  means  of  studying  all  the  modes  simultaneously  is  therefore  to 
isolate  effects  to  a  given  order  of  e. 

For  this  problem  we  need  to  expand  the  wave  function  to  order  e2  so  that 
only  four  time  scales  are  necessary,  i.e.  t-i ,to,ti ,t2.  Of  these  four  time 
scales,  tj  will  make  no  contribution  to  the  calculation  so  it  will  be  eliminated 
from  the  start.  The  time  derivative  in  Eq.  (2.2a)  is  then  written  as 


3  -I  3  ^  3 

3t  3t_j  3t0 


e2 


3 

3t2 


(3.3) 


Once  Egs.  (3.1)  -  (3.3)  are  substituted  in  fiq.  (2.2a)  and  coefficients  of  the 
various  powers  of  e  are  set  equal  to  zero,  the  following  set  of  coupled 
differential  equations  emerge: 


The  initial  conditions  are  then  given  by 

*n(x#t_i  »  0,t0  »  0,t2  *  0)  -  6no  exp(iSx)  (3.9) 

Equation  (3.4)  requires  that  4»0  is  independent  of  -t_1,  i.e.  ^o=^,o(x*to»t2) 
With  this  result  Eq.  (3.5)  can  be  integrated  over  t-i  to  give  an  expression  for 
♦i  in  terms  of  ♦© 

*1  «  iKL^t.j  -  q  cos2x  sin2t.  (3.10; 

In  general,  this  equation  would  contain  an  integration  constant  that  is  a 
function  of  x,  to  and  t2.  Because  this  function  can  be  shown  to  be  a  part  of 
♦o»  it  will  not  appear  in  this  discussion.  The  secular  term  in  Eq.  (3.10)  is 
removed  by  requiring  that  the  coefficient  of  t-i  be  set  equal  to  zero.  This 
condition  yields 

U>°  ~  29  «>s2x  +  i  *  0  (3.111 

and  leaves 

♦l  ■  -iq  cos2x  sin2t_n|)0  (3.12] 

Equation  (3.11)  is  the  one  we  would  have  obtained  if  we  had  followed  the 
traditional  approach  and  averaged  bq.  (2.2a)  over  the  rapid  oscillations. 
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The  solution  o£  Eg.  (3.11)  which  satisfies  the  initial  condition  (3.9)  is 


*0(x,t0,t2)  -  l  o„(t2)me  (x,q)  exp{-ii  (q)t) 


(3.13) 


a  (0)  -  Cp^n  (q) 
n  -2n 


(3.14) 


Here  me  (z,q)  is  a  Mathieu  function  [14]  of  order  a  and  \  (q)  is  its  eigenvalue, 
a  a 

The  constants  C®  (q)  Me  the  Fourier  coefficients  of  me  (z,q),  i.e. 


me  (zfq)  *  l  C®  (q)  exp{i(o  +  2r)z} 


The  substitution  of  Eg.  (3.12)  into  Eg.  (3.6)  and  the  subsequent 
integration  over  t_i  yields  the  following  expression  for 


(3.15) 


,t  ,t  )  =  2q(cos2t  -1)  cos2x^0  ♦  sin2x 
-10  2  -1 


q2 

+  -4-  cos  2x(cos4t  -1)*0  +  B(x,t0,t2)  , 


where  fl  is  an  undetermined  function  that  satisfies  the  initial  condition 


B(x,0,0)  =  0. 


(3.16) 


(3.17) 


There  are  no  secular  terms  this  time.  In  a  similar  manner  Eg.  (3.7)  can  be 
integrated  over  t  after  the  substitution  of  Eq.  (3.16).  Because  we  are  solving 


.'-.■■'.'•I-'---'-- 


I 


for  the  wave  function  to  0(e2)#  we  only  need  to  investigate  the  differential 
equation# 


fa*  a  1  a*0  av0 

- 2q  cos2x  ♦  i  ■*— -  B  =  -i  -r—~  *  8q  cos2xt|»o  +  2  sin2x  -r — 

3z2  «0  3t2  3x 


-  cos2x  -  ♦  3q2(l  -  cos4x)ijio 

3x2 


(3.18 


that  arises  from  setting  the  coefficient  of  t  to  zero.  The  left-hand  side  of 
Eq.  (3.18)  contains  the  same  operator  that  occurs  in  Eq.  (3.11).  We  can 


therefore  write  B  as 


B  -  l  exp{-iX6+2m(q)tT0}. 


(3.19 


When  Eqs.  (3.13)  and  (3.19)  are  substituted  in  Bq.  (3.18)  and  the  coefficients 
of  exp  iXg+2m W  t0  are  set  equal  to  zero#  the  following  equation  arises 

a2  3otn  a 

- 2q  cos2x  +  _  b  =  -i  +  8a  q(cos2x  +  2  sin2x  *— 

3x2  m  B+2nJ  n  3t2 


+  cos2xXg+n)  -  q2(5  +  11  cos4x)  me^ 


(3.2C 


We  can  now  apply  Fredholm's  alternative  theorem  [15]  to  this  equation  which 
requires  that 


(3.23 


6  *  8q[(l  ♦  Xg+2n)  cos2x  ♦  2  sin2x  -  q2(5  +  11  cos4x) 

Equation  (3.21)  with  initial  condition  (3.14)  has  the  obvious  solution. 

°n  *  c-2n"(q)  exP(”iYnC<l)t2)  (3.24 

Equations  (3.12),  (3.13),  and  (3.24)  can  be  combined  by  using  Eqs.  (3.1) 
and  (3.2)  to  give  an  approximate  expression  for  the  wavefunction 

iKx.t)  *  [1  -  icq  sin(2-r/e)]  J  C_2*nCq)®e6+2n(q,x) 

n*  -<*> 


x  exp{-i(Xg+2n(q)  ♦  e2Yn(q))T}. 


(3.25; 


This  wave  function  is  still  properly  normalized  to  terms  of  0(e2).  For  the 
condition  that  we  assumed  at  the  beginning  of  this  section  (e1  »  q  and 
e”1  »  |0|)  the  corrections  to  the  lowest  order  wavefunction  are  negligible.  Ai 
an  exanple  consider  the  terms  in  the  square  brackets.  If  we  calculate  the 
porobability  (2.7)  that  the  electron  will  be  in  a  particular  momentum  state 
after  a  time  t  then  the  second  term  in  the  square  bracket  makes  a  contribution 
that  is  a  factor  (eq)2  smaller  than  the  first  and  can  therefore  be  safely 
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ignored  as  long  as  e  1  »  q.  Even  if  there  were  some  means  of  detecting  this 
small  amplitude  correction ,  we  would  have  to  contend  with  the  fact  that  this 
term  is  oscillating  at  the  frequency  of  the  standing  wave  field. 

The  second  term  in  the  argunent  of  the  exponential  in  bq.  (3.25)  also  does, 
not  appreciably  influence  the  scattering  probability  as  long  as  we  are 
considering  reasonable  interaction  times  (t  <  10).  This  term  introduces  a  slow 
modulation  in  time  to  the  wavef unction.  For  q  «  e”1,  the  period  of  oscillation 
is  much  greater  than  conceivable  interaction  times  so  again  the  correction 
induced  by  the  oscillating  portion  of  the  potential  can  be  neglected.  We  are 
thus  led  to  the  conclusion  that  as  long  as  q  «  e"1  and  | B|«  t  l  the  temporally 
oscillating  term  in  the  vector  potential  can  be  safely  neglected.  This  result 
contradicts  the  conclusions  of  Ref.  [6] . 

The  perturbation  calculation  is  no  longer  valid  when  eq  «  1  since  the 
second  term  in  the  square  brackets  of  Eq.  (3.25)  becomes  ccoparable  with  the 
first  term.  If  we  return  to  the  definition  of  q,  Eq.  (2.9),  and  e,  Eq.  (2.2e), 
then  the  product  can  be  rewritten  as 

2tttoX2I 

eq  *  ~hoi'c  *  2Trr0*2P  <3»26> 

where  I  »  ftwcp  has  been  used  in  the  last  step  and  p  is  the  photon  number 
density.  An  estimate  of  when  the  time  variation  of  the  vector  potential  becomes 
important  is  2nr0X2p  *  1.  In  other  words,  there  is  one  photon  in  the  vicinity 
of  the  electron  at  all  times.  This  condition  makes  sense.  The  time  dependent 
part  of  the  standing  wave  represents  a  stimulated  emission  or  absorption  of  two 
photons  by  an  electron.  This  process  will  not  be  probable  unless  there  are  two 
photons  in  the  vicinity  of  the  electron.  Because  a  photon  cannot  be  localized 
perpendicular  to  direction  of  motion  to  a  distance  less  than  a  wavelength,  the 


vicinity  of  an  electron  is  the  volume  2irroX2.  We  have  therefore  returned  to  the 
condition  given  fay  Gg.  (3.16).  For  X  »  1  urn,  Eq.  (3.16)  requires 
1  «  1014  w/cm  .  Although  these  intensities  can  be  achieved  for  focussed  laser 
beams  they  cannot  be  achieved  in  the  1  cm  beams  needed  to  generate  a  reasonable, 
interaction  time  for  the  electrons.  On  the  other  hand,  C02  lasers  (X  *  10  yon) 
require  1  a  10llW/cn2  forEq.  (3.16)  to  be  satisfied.  This  is  presently  possible. 

When  the  small  terms  in  Eq.  (3.25)  are  dropped,  the  wave  function  assumes 
the  form 

» 

'l'(x.t)  «  l  c!2n^meB+2n(x,<^  exP{"iXg+2nCq)Tj  •  (3.27) 

n«-°» 

The  probability  for  scattering  from  an  initial  state  with  momentum  p^  ■  ehk  to  a 
final  state  with  momentum  =  (8  +  2r)hk  after  a  time  t  is  then  given  by  Eq. 
(2.7)  as 

PrCB,T)  =  |  l  c^n(q)  exp{-iX6+2n(q)T>|2 

n=-°° 

(3.28) 


and  satisfies  Eq.  (2.8).  Although  very  different  in  form,  the  continued 
fraction  expression  for  the  scattering  amplitude  derived  in  Ref.  [5)  is 
identical  to  Eq.  (3.27).  The  advantages  of  the  Mathieu  function  expansion  over 
its  continued  fraction  form  is  that  Eq.  (3.28)  is  easier  to  approximate 
analytically  and  evaluate  numerically.  Because  the  Fourier  coefficients  of  the 
Mathieu  functions  satisfy  a  three  term  recursion  relation  [14],  Mathieu's 
equation  can  be  written  as  a  matrix  eigenvalue  problem  where  the  matrix  is 


tridiagonal.  The  eigenvalues  x.  _  and  the  Fourier  coefficients  c?+2n  are  then 

.  p+zn  zr 


easily  found  numerically  using  any  standard  program  to  diagonalize  the 
tridiagonal  matrix. 


Figure  2  shows  the  scattering  probability  as  a  function  of  q  for  pi  ■  hk; 
p  »  hk,  -hk,  3tik,  -3hk  and  t  ■  2  for  a  standing  wave  with  a  wavelength  of  1 
vm.  The  interaction  time  is  chosen  so  as  to  represent  the  approximate  time  a 
2U0  eV  electron  takes  to  travel  one  centimeter.  These  values  of  the  parameters  . 
are  chosen  as  representative  of  a  feasible  experiment.  One  noticable  feature  of 
these  graphs  is  the  increasing  sensitivity  of  the  scattering  probability  to  the 
standing  wave  intensity.  This  result  is  not  surprising  since  increasing  q 
increases  the  nuriber  of  terms  that  contribute  to  the  summation  in  Eg.  (3.13). 
Each  additional  term  adds  another  and  higher  frequency  component  to  the 
scattering  amplitude.  Evan  for  relatively  low  intensities  (such  as  at  the  peak 
of  the  first  maximum  of  the  probability  for  scattering  from  hk  to  -hk  (Fig. 

(lc) ) ,  an  increase  of  q  by  a  factor  of  2  is  sufficient  to  move  the  scattering 
probability  to  near  zero.  This  sensitivity  could  explain  the  difficulty  with 
seeing  this  effect  in  the  early  experiments  (7-12).  The  variation  of  the 
probability  with  q  becomes  less  severe  if  0  is  increased  and  higher  order 
scattering  is  considered.  This  result  is  shown  in  Fig.  3.  Although  the  peak 
scattering  probability  is  reduced  from  the  case  shown  in  Fig.  2  it  is  still 
respectable. 

The  probability  for  scattering  from  hk  to  -hk  and  from  3hk  to  -3hk  as  a 

function  of  the  interaction  time  t  is  shown  in  Figs.  4  and  5,  respectively,  for 

several  values  of  q.  Again,  an  increase  in  q  causes  the  probability  to 

« 

fluctuate  more  rapidly  but  this  time  as  a  function  of  t. 

The  last  parameter  that  can  be  easily  varied  in  an  experiment  is  the  angle 
between  the  electron  beam  and  the  axis  of  the  standing  wave  or  equivalently  B. 

In  Fig.  6  the  scattering  probability  for  the  transition  from  Bhk  to  (B-2)hk  as  a 
function  of  0  is  shown  for  two  different  values  of  q.  These  graphs  show  clearly 
that  for  this  particular  transition  the  scattering  probability  is  synraetrical 
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about  the  point  8*1.  However,  8  ■  1  is  not  necessarily  the  point  of  maxima 
scattering.  This  could  also  explain  the  failure  of  several  of  the  experiments 
to  observe'Kapitza-Dirac  scattering. 

IV.  STRONG  COUPLING  AND  HIGH  ENERGY. 

As  was  pointed  out  in  Section  III,  the  MTS  analysis  that  gave  the  leading 
behavior  in  the  form  of  slowly  modulated  eigenfunctions  of  the  time-averaged 
problem  breaks  down  if  q  *  0(l/e) .  To  examine  this  case,  we  employ  a  variant  of 
the  WKB  method  that  allows  for  the  possibility  of  resonance,  caused  by  the  time 
dependence  of  the  potential.  We  find  that  if  certain  conditions  are  met  between 
the  wave  characteristics  of  the  incoming  electron  wavefunction  and  the  standing 
wave  field,  the  amplitude  evolves  to  very  large  localized  maxima  near  focal 
points  and  caustics  of  the  rays  of  the  corresponding  classical  system.  We 
expect  that  near  such  caustics  the  effects  of  many-particle  interactions  and 
self-radiation  will  become  important.  A  realistic  analysis  of  these,  including 
relativistic  effects,  will  be  presented  in  a  subsequent  paper. 

This  section  is  organized  as  follows:  In  part  A  we  discuss  the  expansion 
used  and  outline  the  calculation.  In  part  B  we  present  a  perturbative  treatment 
of  the  Hamilton-Jacobi  equation  for  the  rays;  a  lowest  order  resonant  case  is 
investigated  in  detail,  and  found  to  lead  to  the  focusing  of  rays  and  caustic 
formation.  Finally,  in  part  C  we  analyze  the  effects  of  focusing  in  the 
classical  problem  on  the  probability  amplitude. 

A.  The  Quasiclassical  Expansion 

Although  problems  where  a  high  frequency  appoximation  is  relevant  have 
been  studied  extensively  for  time-independent  potentials  very  few  results  of 
this  type  exist  for  the  non-separable  time  dependent  case  [16] .  In  our 
discussion  of  the  Kapitza-Dirac  problem  we  shall  use  the  potential 
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V(X,t)  «  COS2  X  cos2 t 

but  the  method  can  be  used  for  more  general  time-periodic  potentials.  For 
problems  of  this  type  we  expect  quasiperiodic  behavior. 


In  the  previous  section  we  saw  this  quasiperiodicity  arise  in  our  MIS 
treatment  which  was  valid  (at  least  formally)  for  moderate  energies.  In 
examining  the  high  energy/high  frequency  behavior  in  a  system  without  internal 
degrees  of  freedom,  whose  classical  counterpart  is  described  by  a  Hamiltonian 
H(p,q,t)  the  ansatz 

♦  *  A  exp(i S/h)  (4.1 


is  used.  Substitution  in  the  Schrodinger  equation 


a<ji 


H(p»q#t)  i|>(q,t)  «  -  i  ~ 


(4.2 


(where  p  =  -  j~-  •—)  results,  to  leading  order  in  h_1  in  the  equation 


(4.3) 


i.e.  the  classical  Hamilton-Jacobi  equation,  showing  that  the  phase  S 
corresponds  to  the  classical  action.  Ihe  next  order  produces  an  equation  for 
the  amplitude  (*  «  A2): 


♦ 


l  TZ  (v*)  -  0 


(4 


whose  characteristics  (rays)  are  the  same  as  those  of  (4.3).  From  (4.4)  it  is 
seal  that  it  is  concerved  in  ray  tubes  so  that  if  a  tube  collapses  ir  (and  A) 
becomes  infinite.  This  happens  on  caustics  of  the  Hamilton-Jacobi  equation 
which  are  envelopes  of  families  of  rays  in  the  (q,t)  plane.  If  we  think,  in 
(q,t,Sq)  space,  of  the  surface  formed  by  the  rays  through  sane  initial  curve 
(e.g.  if  S^(q,0)  is  given),  than  the  caustics  are  the  singularities  of  its 
projection  on  the  (q,t)  plane,  corresponding  to  folds,  etc.  of  the  surface. 

It  must  be  understood,  as  was  shown  by  Buchal  and  Keller  [17] ,  that  the 
higher  order  terms  neglected  in  (4.4)  will  becane  large  at  the  caustic  and  thus 
they  must  be  included  there.  This  effectively  gives  rise  to  a  boundary  layer  in 
the  vicinity  of  the  caustic  in  which  the  amplitude  is  large  but  still  finite. 

The  thory  of  geometrical  optics  allows  us  to  connect  through  a  caustic  by 
including  an  appropriate  phase  shift  in  the  (coriplex)  amplitude.  The  value  of 
this  phase  shift  is  m(v/2)  where  m  is  the  order  of  degeneracy  of  the 
projection,  or  equivalently,  the  nunber  of  dimensions  lost  by  the  ray  tube  at 
the  caustic  or  the  change  in  the  number  of  branches  in  the  vicinity  of  the  point 
considered.  The  sign  is  chosen  according  to  whether  the  caustic  is  traversed  in 
the  direction  of  increasing  (-)  or  decreasing  (+)  S  [18] .  This  was  first 
realized  by  J.  Keller  in  his  1958  paper  [18]  where  he  also  pointed  out  the  need 
for  many-branch  descriptions  of  the  form 

r 

*  *  l  \  exptiS^/h) 


(4.5) 


where  r  is  the  number  of  rays  through  the  point  under  consideration.  Demanding 
single  valuedness  for  <J»  he  showed  that  in  nonseparable  systems,  quantum  nanbers 
for  bound  states  are  in  general  quarter  intergers  (asymptotically) .  V.  Maslov 
[19]  proved  the  asymptotic  character  of  this  approximation  as  h  -*■  0  for  a 
special  class  of  problems:  essentially  "nice"  time-independent  potentials  in 
several  space  dimensions  and  initial  conditions  that  vanish  outside  some  finite 
region.  He  showed  that  for  ESq.  (4.2)  with  H  in  the  form 

H(p,q)=  j  l  p2  +  V(q)  (4.6) 

and  with  initial  conditions  of  the  form 


i|i(q#0)  *  a(q)  exp(ib(q)/h) 


(4.7) 


with  a(q)  zero  outside  a  finite  region  that  an  asymptotic  expression  for  \l>  is 
given  by 


i|»(q#t) 


r 

l 

K=1 


a(qk)J(q?qok) 


-h 


exp 


(^sk(,,t>  -  i|vk)  ♦  Olh) 


(4.8) 


here  J  is  the  Jacobian  of  the  mapping  from  the  initial  point  q0  of  the  k*th  ray 
through  q  to  the  point  q  induced  by  the  classical  Hamiltonian  flow,  S  is  the 
classical  action  along  that  ray  and  is  the  Morse  index  of  the  ray,  equal  to 
the  sum  of  the  numbers  m  for  every  encounter  of  the  given  ray  with  caustics. 
This  expression  fails  at  caustics,  where  the  Jacobian  becomes  singular. 

In  part  B,  we  derive  a  formula  similar  to  (4.8)  for  our  problem,  without 


attempting  a  rigorous  justification  (for  which  we  would  have  to  replace  our 
plane  wave  initial  conditions  with  a  function  equal  to  it  inside  some  finite 


region,  that  is  large  compared  to  the  wavelength,  and  zero  outside) .  The 
classical  problem  for  the  Kapitza-Dirac  potential  corresponds  to  a  pendulum  with 
oscillating1  strength,  described  by  the  Hami 1 ton-Jacobi  equation 

St  +  S2  +  e4q  cos2  x  cos2t  ■  0  (4.9) 

with  initial  condition 

S(x,0)  -  e26x,  (4.10) 

Here  6  is  the  wavenumber  of  the  incoming  plane  wave.  We  consider  especially 
8*0  (1/e2),  q  9  0(l/e2).  It  can  be  seen  that  when  6  is  in  a  rational  relation 
with  the  periods  of  the  potential  the  classical  problem  exhibits  resonance  (at 
least,  in  the  sense  of  perturbation  theory),  resulting  in  strong  bending  of  the 
rays.  We  treat  the  case  B  •  1/2 e 2  (lowest  order  resonance)  which  leads  to 
caustic  formation  in  a  time  scale  of  order  1/e. 

The  leading  behavior  of  the  rays  if  found  as 

x  »  t  ♦  Co(t,C)  +  0(e) 

P  -  i-e(icoT]  +  0(e2)  (4.11) 

where  C0  is  a  solution  of  the  pendulum  equation 
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2Cqtt  +  sin2Co  ■  0 


C0 (0,5)  -  C#CoT (0,5)  -  0 


(4.12) 


and  £,  fixed  along  a  cay,  gives  its  position  at  t  ■  0  (t  »  et  is  a  slew  tine 
scale).  We  see  that  the  rays  ace  split  into  groups  performing  an  oscillation 
around  centers  which  drift  with  uniform  velocity.  It  is  clear  from  (4.11)  and 
(4.12)  that  caustics  (locations  of  crossing  rays  characterized  by  the  condition 
x  *  0)  are  cusp-like.  They  form  at  the  center  of  each  group.  (Fig.  7)  and  move 
out,  toward  the  saddle  points.  New  cusps  form  periodically  (with  period 
T  “  2k/ t).  Every  cusp  narks  an  S-shaped  folding  of  the  surface  P  ■  P(x,t)  (Fig. 
8)  (caustics  are  singularities  of  the  projection  of  this  surface  on  the  x-t 
plane,  i.e.  they  correspond  to  Px  becoming  singular).  The  amplitude  is  found  to 
leading  order  as 


A0(T,5>  -  (Coe)"^  +  0(e) 


(4.13) 


On  caustics  x^  »  0  implies  Cog  *  0  and  A  appears  to  become  singular.  The  reason 
is  easy  to  see  if  we  rewrite  the  amplitude  equation 


*  e/  *  ic\x  ‘  0 


A(x,0)  -  1 


in  characteristic  coordinates,  (t,c),  in  which  it  becomes 


(4.14) 
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felt  I  A-° 


(4.15) 


Clearly,  the  0(e2)  tern  which  was  neglected  in  deriving  (4.8)  becomes  important 
near  the  caustic.  Ibis  suggests  studying  Eg.  (4.15)  by  techniques  of  singular 
perturbation  theory.  This  approach  was  taken  for  the  reduced  wave  equation  by 
Buchal  and  Keller  117].  It  produces  the  correct  leading  behavior  near  a  smooth 
convex  caustic. 

An  expansion  for  <j>  is  utilized  in  the  form 


*  *  l  ^  expUS^/e2) 


(4.16) 


with  k  ranging  over  the  number  of  distinct  branches  of  the  surface  P  *  P (x,t) . 

Brandies  join  at  caustics.  A  local  coordinate  system  is  introduced  at  caustics 
and  Eq.  (4.15)  is  scaled  to  give  the  leading  behavior  in  their  vicinity. 

Away  from  cusps,  two  branches  join  at  a  caustic,  an  incoming  and  an 
outgoing  branch.  Using  a  matching  procedure,  we  can  determine  A  on  the  outgoing 
brand)  from  the  (presumed  known)  values  of  A  on  the  incoming  branch.  To  carry 
out  the  matching  it  is  necessary  to  analytically  continue  A  into  the 
(classically  forbidden)  side  of  the  caustic  inaccessible  to  the  rays  and  demand 
that  the  amplitude  decreases  exponentially  away  from  the  caustic  in  this  region  (Fig.  9) 
The  main  result  of  this  analysis,  applied  to  our  problem,  is  that  near  the 
caustic  the  amplitude  is  of  order 


A  %  0(c"1/6) 


(4.17) 
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which  night  play  a  very  significant  role  in  situations  where  there  exists  the 
possibility  of  nonlinear  interaction  (e.g.  refractive  media,  radiation,  etc.). 
The  singular  perturbation  analysis  of  Eq.  (4.15)  can  be  carried  out  for  other 
cases,  e.  g.  near  the  "cusp"  (which  only  appears  as  a  cusp  in  the  x-t  plane,  being 
a  smooth  curve  on  the  surface  P  -  P(x,t)  itself)  or  even  for  higher  dimensional 
caustics.  One  can  use  this  approach  to  get  leading  order  estimates  for  the 
amplitude.  To  carry  out  the  calculation,  especially  in  higher  dimensional 
problems,  one  can  be  guided  into  the  proper  scaling  by  considering  the  normal 
forms  that  exist  giving  generic  coordinate  systems  of  the  surface  near  various 
types  of  caustics  (for  a  classification  and  a  list  of  normal  forms  see,  e.g. 
Arnold  [20]).  This  way  we  can  find,  for  example,  that  near  the  cusp 

A  -v  0(e'1/4)  (4.18) 

B.  The  Expansion  Away  From  Caustics 

In  Sec.  Ill  we  were  led  to  study  a  problem  of  the  type 

-  qV(x,t'/e2)  -  -  ii|» 

**  £  f 


i|*(x,0)  »  exp(iex) 


Rewritting  in  terms  of  t  *  t 'A  »  Bq.  (4.19)  becomes 


e2(^xx  -  <jV(x,t)^)  «  -  ii|»t 
(x,0)  *  exp(iBx) 

To  keep  the  discussion  general ,  we  leave  q  as  a  free  parameter. 
6  has  the  form 

6  “  ~z  (6q  +  a®l  +  °2&2  +  •••) 


(4.20; 


We  assure  that 


(4.21] 


where 


a  - 

In  Sec.  Ill,  the  case  q  *  0(1),  6  ■  0(1)  was  studied.  Here  we  extend  the 
discussion  to  other  interesting  limits.  Using  the  ansatz 

~  I  \  exp(iSk/c2)  (4.22] 

k 

we  find  that  each  branch  must  satisfy  (for  simplicity  we  drop  subscripts) 


-  tL.  (S2  ♦  S  +  e4qV)  -  i  (A  +  2S  A  +  S  A)  +  e2A  =  0 


(4.23 


*  •  * 


We  recognize  the  expression  in  the  first  parenthesis  as  the  Hamilton-Jacobi 
equation  for  the  classical  action  and  we  resolve  the  apparent  ambiguity  in  A,S 
by  requiring  that  S  is  precisely  the  classical  action,  so  that 


St  +  Sx  +  c4<^(x#t)  *  0 


(4.24) 


Then  the  amplitude  A  satisfies 


A.  +  2S  A  +  S  A  -  ie2A  «  0 
t  x  x  xx  xx 


(4.25) 


with  initial  conditions 


S(x,0)  *  e28x,  A(x,0)  *  1 


(4.26) 


It  is  easily  seen  that  the  subcharacteristics  of  Eq.  (4.25)  found  by  neglecting 
the  term  -ie2Axx  are  the  same  as  the  characteristics  (rays)  of  Eq.  (4.24). 

Since  the  phase  equation  is  independent  of  the  amplitude  it  is  studied  first,  to 
determine  the  ray  structure  of  the  problem. 

Utilizing  t  as  the  parameter  along  the  rays  and  £  as  the  ray  variable,  i.e. 
assuming  that  the  equation  of  a  ray  is  x  *  x(t,£)  with  x(0,£)  «  £,  we  rewrite 
Eq.  (4.24)  in  characteristic  form 


*  2P 


(4.27a) 
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(4.27b) 


*  2P2  +  Q  «  P2  -  o2V 


(4.27c) 
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|fi--  C2V 
at  t 


(4.27d) 


where  we  let  *  P,  S{  8  I)  and 


3_  _  3_  3_  3x  3_ 

at  =  at  “  at  x  at  e  ax|t 


is  the  directional  derivative  along  the  rays.  The  x-P  system  (Eg.  (4.27a,b)  can 
be  solved  first ,  then  (Eq.  (4.27c)  and  Eg.  (4.27d)  to  allow  us  to  find  S. 

He  shall  study  Eg.  (4.27a,b)  using  MTS  perturbation  theory.  Ms  expect  the 
weak  (0  (a2) )  forcing  to  result  in  a  slow  modulation  of  the  rays.  The  time  scale 
for  the  lowest-order  resonant  case  we  shall  consider  is 


crt,  o  88  e2q^ 


(4.28) 


while 


X  *  x0(t,T#£)  +  oX1  +  o2x2  ♦  0(o3) 


P  «  P0  +  oPx  +  02P2  +  0(o3) 


with  initial  conditions 


x(0,£)  -  £ 


P<0,£)  -  e26  *  60  *  a&0  +  °2H  + 


*_*  •  *  •  •_ a\ji  "j*  *  •  “  * 


(4.29) 


(4.30) 


>  J.  -  .  •  •  •*  ■  .  V.  -  .  *  »%  •  *  .  •  «  v.  '  •  •  .  *  .  •  .  “  -  ‘  *  >>  •  *v.  v*\*  \ 
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Substituting  the  expansions  (4.29)  into  (4.27)  and  equating  coefficients  of  like 
powers  of  e  we  find  (recall  that  a/at  ■  a/at  ♦  eU/dt))  the  hierarchy  of 
equations 


0(1)  Xot  -  2Pfl 


P  «  0 
ot 


(4.31) 


0(o)  xu  «  2Pj  -  xQt 


P  -  -  P 

It  OT 


(4.32) 


0(o2)  x  «  2P  -  x 

2t  2  OT 


with  I.  C. 


P  -  -  P  -  Vx (x  ,t) 

2t  It  *  0 


(4.33) 


x0(0,£)  *  *i(0#O  ■  Of  i  ■  1,2,  ... 


Pi(0,t)  ■  6i#  i  *  0,1,2, 


Solving  we  find 


0(1)  Xq  *  2gjt  +  Cq (t ,£) ,  C0(O,e)  ■  £ 


(4.34) 


Pq  *  Bq(t,C)»  ®0^0,C)  ■  Bo 


(4.35) 
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(assuming  the  limit  exists,  which  it  does  for  V  periodic  in  t).  The 
construction  can  be  continued  to  higher  order,  but  here  we  only  need  the  leading 
behavior.  We  found 

x  -  230t  +  C0(t,5)  +  0(o) 

P  -  Bo  +  otjC0T)  +  0(o2)  (4.40) 

where  Co  satisfies  Eq.  (4.39).  For  the  potential 

V(x,t)  *  cos2x  cos2t  -*•  V  *  -  cos2t  sin2x  (4.41) 

we  have 

|  Vxdt  *  "  88o  C°S2(2B°t  +  Co)  "  8(2T0±lT  COS2l(2e±1)t  +  Colf 
o  0 

a  bounded  function  if  So  t  0,  ±  1/2.  In  this  case,  Bq.  (4.39)  gives 
C0U,t)  *  Sit  +  C  and  the  slow  time  scale  t  only  appears  as  a  modifications  of 
the  fast  scale  t.  In  general  for  a  problem  of  this  kind,  we  should  consider  a 
modified  fast  time  scale 

t*  m  w(ct)  t  *  (o)g  +  0<J)J  +  »..)t 

and  a  slow  time  scale 

t  *  o°t 


Where  a  is  determined  by  the  type  of  resonanoe  expected.  To  treat  the  simplest 
case,  we  consider 


^0  *  2  '  *  0  i  *  1,2,  ... 

q  *  e~2,  a  ■  e  (4.42] 

Now  we  have  a  lowest  order  reronance;  the  integral  in  Eq.  (4.39)  becomes 

J  Vxdt  -  i  cos(2t  +  2C0)  +  cos(4t  +  2C0) 

0 

-  cos2C0  -  J  sin2C0  (4.43) 

and  the  equation  reads 

2CQtt  +  sin2C0  =  0  (4.44) 

with  l.C. 

C0(O,5)  «  Cot(0,C)  =  0  (4.45] 

Finally,  utilizing  Eq.  (4.40)  in  Eq.  (4.27c)  we  find 

s“it  +  ic°  +  °(e) 

Turning  our  attention  to  the  amplitude  equation 


V.V.V  *•*. 


■ 


A.  +  2PA  +  P  A  -  c2iA  -  0 

t  X  X  XX 


tie  transform  it  to  characteristic  coordinates 

3  -1  3  3  __  a  3 

3x  ‘  x?  ’  3t  2P  3x  ’  35  ( 

to  find 

\  *  *;v  - e2i  [xv  kY » •  °  < 

Assuming 

A  *  A0  +  eAx  +  0(e2) 

and  using 

-  c,t  ♦  «(«» 

p{  ■  £<Kte>  *  °<‘2> 

we  have 

\  *  (Coc  +  0(c))  (e(icOTC)  +  0(e2))A  -  0(e2)  (- 

from  which  we  find 


-  .  V  ’  •••  •• 
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0(1)  A  «  0  •*  A0  «  Ac(t#c)  with  A0(O,5)  «  1 


°<‘>  *»  *  V  +  £  ■  1  Cou  ■  0 

ot 


A,  *  -  (A  +  I  A  +  (bounded  part) 

1  Ot  2  C 


OC 


so  finally,  demanding  that  A^be  bounded  we  find  that  A0  is  given  by 


A0(t,£)  -  (Co5) 


(4.50) 


Collecting  our  results  we  see  that  the  leading  asymptotic  behavior  of  is  given 


as: 


♦  ^  ((C„r)",S  +  0(e))  exp  -i-  (lt  +  C0(tK)  +  0(e))  (4.51) 

o?  2e2  2 

with  C0  satisfying  Eq.  (4.44)-(4.55) .  We  see  (Pig.  7)  that  the  resonance  causes 
a  slow-time  bending  of  the  rays  so  that  the  surface  P  *  P(x,t)  given  by 

x  «  |t  +  C„(t,0  ♦  0(e) 

P  «  c(jC0t)  +  0(e2)  (4.52) 

develops  folds  at  points  where  X  ■  0,  or  equivalently,  x  *  0  -*■  C  *  0. 

P  5  Ot, 

Clearly  this  vanishing  leads  to  a  blowup  in  Eq.  (4.50) ,  invalidating  our 
asymptotic  solution  at  such  points.  What  went  wrong  in  our  analysis  is 
apparent:  Eq.  (4.48)  was  analyzed  by  a  regular  perturbation  expansion. 


,\-Vv 

.V. 


•*“  m  •  ,*»  ,**  ,*»  *»  *• 

V  •  «  «  .*,«,•  ,  •  .  r< k  .  ,  ,  «•  . 
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is  • 


affectively  neglecting  the  0(e2)  term  in  finding  A  .  In  the  neighborhood  of 
envelopes  of  rays,  where  x^  -  0,  this  teem  clearly  will  be  important  so  that  A 
would  stay  finite,  although  it  will  become  very  large  there.  For  this  reason 
points  where  x^  »  0  are  called  caustic  points.  In  the  next  section  we  complete 
the  picture  that  this  analysis  presents  by  discussing  the  behavior  of  solutions 
in  the  neighborhood  of  caustics. 

C.  The  Behavior  Near  Caustics 

We  now  investigate  the  behavior  of  the  amplitude  near  caustics. 
Following  Ludwig  [21],  we  recognize  expansions  of  the  form 


*  *  I  \  exp(iSk/e2) 


as  resulting  when  expanding  integrals  of  the  type 


(4.53 


J  A(x,t,z)  exp(iB(x,t,z)/e2)dz 


by  the  method  of  stationary  phase.  Clearly,  assuming  that  can  be  represented 
as  a  superposition  of  plane  waves,  we  will  have,  as  e  -*■  0,  contributions  from 
points  where 


B  «  0,  B  f  0 
z  zz 


(4.54) 


Such  points  can  be  shown  to  correspond  to  regions  away  from  caustics.  Where  the 
expansion  (4.8)  is  valid.  Smooth  caustics  correspond  to 


B  -  0,  B  4  0 
zz  zzz 


(4.55) 


y-  .y-y-.- .*• '.-.V-'I-.V-V  Sv-V-'V.-V-V. 


while,  finally  cusped  caustics  are  found  if 


B 

zzz 


0, 


B 

zzzz 


f  0 


(4.56) 


Systems  of  rays  in  two  dimensions  have  in  general  singularities  of  the  "fold" 
(i.e.  smooth  caustics)  or  "tuc"  (i.e.  cusped  caustics)  kinds  (see  Arnold  [20]), 
unless  same  exceptional  situations  happen.  An  easy  calculation  shows  that  for 
our  problem  the  worse  thing  that  can  happen  is  a  cusp,  so  we  only  need  to 
consider  the  cases  for  Egs.  (4.55)  -  (4.56) .  He  can  encompass  all  cases  by  the 
exact  change  of  variables 


B(x,t,z)  -  S(x,t)  +  ri£  -  r2C2  ♦ 


(4.57) 


(with  ri  *  r^(x,t))  which  is  related  to  the  normal  form  that  obtains  near  a 
"tuc"  singularity  (Gq.  56)).  Then,  by  following  the  same  procedure  as  in  Ref. 

[21]  for  the  Schrodinger  equation,  we  can  construct  uniformly  valid  expansions 
which,  away  from  caustics,  reduce  to  Eq.  (4.8),  near  smooth  caustics  have 
leading  terms  involving  the  Airy  function  and  its  1st  derivative  and  near  the  cusp 
involve  a  generalized  Airy  function  with  two  arguments  and  its  two  1st 
derivatives.  By  considering  the  asymptotic  behavior  of  the  Airy  functions  we 
can  find  the  order  of  the  amplitude  in  the  vicinity  of  caustics,  and  by 
analytically  continuing  to  the  classically  inaccessible  regions  we  see  that  our 
solutions  decay  exponentially  there.  Since  here  we  are  not  as  much  interested 
in  the  detailed  behavior  at  all  regimes  but  mainly  on  the  order  of  magnitude  of 
the  focusing  at  caustics  we  shall  follow  instead  a  singular  perturbation 
approach  due  to  Buchal  and  Keller  [17].  He  introduce  boundary  layers  at 
caustics  and  determine  the  solution  in  their  vicinity  by  stretching  coordinates 
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and  matching  with  the  "outer"  solution.  We  shall  use  this  method  to  give  the 
expansion  near  a  smooth  convex  portion  of  the  caustic  and  then  we  shall  give  an 
argument  about  the  behavior  near  the  cusp.  We  prefer  the  singular  perturbation 
approach  because  it  would  be  applicable  in  a  more  interesting  setting,  namely  a 
weakly  nonlinear  equation  that  has  singularities  related  to  folding 
subcharacteristics.  Clearly,  in  such  cases,  it  would  be  impossible  to  argue 
from  the  point  of  view  of  integral  representations,  while  the  boundary  layer 
idea  would  still  be  relevant  (after  suitable  modification) . 

We  wish  to  study  the  equation 


PF 

A^  +  -2-  A  -  ie2 


J_  d_ 

h  \ 


(4.58) 


near  characteristics  where  =  0  and  the  second  derivative  term  becomes 
singular  and  cannot  be  ignored.  The  equation  of  a  ray  was  found  to  be  (to  the 
leading  order) 


x  -  t  +  C0(£,t)  +  0(e)  (4.59) 

while  on  caustics  x^  *  0,  i.e.  CQ^  *  0.  We  let 

C  _  (C,t)  *  0  -*■  T  *  b(£),  i.e.  t  *  I  b(£) 

Ot.  c 

Then  the  equation  of  a  caustic  will  be 
x  -  I  b(5)  +  C0(£,b(£))  +  0(e) 


(4.61) 


With  Ref.  [17]  we  introduce  a  new  coordinate  system  (r,s)  where  r  is  arclength 
along  the  caustic  and  s  is  arclength  along  the  ray  tangent  at  r,  measured 
respectively,  from  sane  fixed  point  and  from  the  point  of  contact  (positive 
after,  negative  before) (Fig.  10) .  An  easy  calculation  shows  that,  near  the 
caustic  (i.e.  for  t  -  b(£)  "small") 


s 


-  b(O)  +  oa) 


(4.62) 


Along  the  caustic  we  find  (for  Co  sane  fixed  number) 


(b(£)  -  b(C0>)  +  0(1) 


(4.63) 


while  the  radius  of  curvature  of  the  caustic  p  is  found  from 


i  e2<coct  +  ¥W 

p  b5[l  *  (1  ♦  cC0t)2]3/2 


(6,b(C)) 


(4.64) 


where  C0^T  <  0,  b  >  0  on  the  caustic  (Fig.  11).  The  analysis  we  present 
here  will  be  valid  unless  p  ■  »,  which  happens  at  the  inflection  point.  Tfte 
difficulty  is  not  a  real  one  though:  it  can  be  resolved  by  using  a  better 
coordinate  system,  e.g.  arclength  along  the  caustic  and  normal  distance  from  the 
caustic  (our  present  coordinate  system  would  become  rather  awkward) .  In  the 
convex  part  of  the  caustic,  we  set 


p«—  RG),  (R  *  0(1)) 
e2 


(4.65) 


Utilizing  Bgs.  (4.62)  -  (4.63),  we  find  that  the  equation  becomes  (to  leading 


<£a  ♦  S-K.  -LS?_L|_fiAl 

S  -  /2s  C2  s  ^  3s  U  sj 


(4.66) 


all  the  terms  of  this  equation  will  be  of  the  same  order  if  we  introduce  a 
scaled  variable 


o  -  e2/3s 


Clearly  (see  Eq.  (4.62),  o  is  another  slow  time  scale,  o  e2/,3(t  -  bCcD/e ) .  In 


terms  of  it  we  have 


2A  .0 

a  a  o  3a  (o  J 


(4.67) 


This  equation  is  analyzed  in  Ref.  [17]  where  a  matching  with  the  outer  expansion 
is  carried  out  to  leading  order.  Adapting  the  results  to  our  case,  we  have  that 
the  outer  solution  is 


*o  *  |Co5f%  %  "  b<C>)C0tT(S,b)|',s,  T  <  b(C) 


A„  -  exp(-iri/2)  |C0£|  ,  t  >  b(C) 


(in  agreement  with  Eq.  (4.8),  and  the  discussion  in  part  A),  while  the  inner 
solution  is 


*  AU-(S)2/3o2R'2  «xp<i  2li)l-«xp(-i  — ) 


37 


(according  to  t  -  b  J  0) . 

Therefore  we  can  say  that  the  caustic  boundary  layer  is  of  thickness  e“'  "’(in  t) 
and  the  amplitude  there  is 


.2/3.. 


A  -  0(e'1/6) 


We  are  able  to  carry  the  matching  out  between  incoming  and  outgoing  branches  by 
superimposing 


l b  *  \b.  ♦  \/) 

T  rincommg  Toutgomg 


and  demanding  that  the  continuation  in  the  classically  inaccessible  region  is 
exponentially  decaying  away  from  the  caustic. 

It  must  be  noted  that  we  did  not  give  a  uniformly  valid  expansion  but  only 
leading  terms  near  and  away  from  the  caustic.  As  shown  in  Ref.  [21],  a 
uniformly  valid  expansion  contains  a  term  that  in  the  caustic  is  of 
involving  the  derivative  of  the  Airy  function,  which  in  some  intermediate  region 
becomes  important  and  must  be  included.  A  similar  term  should  be  obtainable  by 
the  matched  expansion  method  if  we  carry  it  to  the  next  order. 

we  close  this  discussion  by  giving  an  argument  for  the  order  of  A  at  the 
cusp.  By  introducing 


T  “T0  , 


So 


where  (tq,£o)  is  the  position  of  the  cusp  with 
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*0  "  ~  +  CoUo/Tq) 


we  have  that  in  the  neighborhood  of  the  cusp  rays  are  given  by  (to  leading 
order) 


X  -  xn  *  I  +  C  ZT  +  Ic  Z3 
0  *  e  0£t  6  0£££ 


(4.68) 


(all  functions  evaluated  at  (to»£q) ) •  Then 


*(*  V  +  Kwz2  +  - 


VeTC.«  +  - 


It  can  be  easily  seen  that  at  (50,t0)#  (C0£t<  °»  Co5C5  >  °»  50  that  Eq.  (4.68) 

is  a  surface  in  (X/Z,T)  space  that  forms  a  fold  at  T  *  0  as  expected. 

\ 

Transforming  the  amplitude  Eq.  (4.58)  to  (Z,T)  we  have  (letting  C  *  a, 

C  *  b) 


rA_  +  —  ea  A  -  ic2  1  3 

*  2  aT  ♦  hbZ2  aT  ♦  HbZ2  32 


A-0 


and  all  terms  are  of  the  same  order  if  we  introduce  scaled  variables 


V  •  e2°a,  Z  ■  e°c 


where 


«*  •*,  /  •  ‘ ,  •*.  •*.  •  ,  **,•'.*  •  •*.  *  \  •*.  *.  •*. 
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o  ■  1/4 

He  note  that  a o  +  Jsbc2  -  0  giving  the  caustic  is  a  smooth  curve:  its  projection 

in  the  x-t  plane  has  a  cusp  at  a  ■  ;  ■  0  (x  ■  x0).  The  outer  expansion 

A«  *  (C  )~i/2  becomes,  in  the  vicinity  of  the  cusp 
o  05 

Ac  *  <C05tT  +  J  z2) 

and  in  terms  of  inner  variables 

A0  *  e'1/4|aa  +  \  b?2|‘1/2  (4.69) 

% 

Without  being  proof  that  A  «  0(e-1^4)  near  the  cusp,  it  still  provides  us  with 
a  good  indication  as  to  what  we  can  expect  there  (it  turns  out  to  be  what  we 
would  get  if  we  carried  out  Ludwig's  method  for  this  case)  [22]. 

V.  CONCLUSION 

In  this  paper  we  have  investigated  Kapitza-Dirac  scattering  in  two 
experimentally  interesting  regimes  of  values  of  the  governing  parameters  q,  g 
and  e.  The  choice  of  two  regimes  was  based  on  the  distinct  methods  used  for 
their  study  as  well  as  qualitative  differences  in  the  behavior  of  the  system. 
Previous  analyses  of  this  problem  let  to  conflicting  answers  on  the  role  of  the 
time  dependent  portion  of  the  Hamiltonian.  In  the  first  case  examined  in  this 
paper  (q,  g  «  e"1)  it  was  found  that,  to  lowest  order,  treatment  of  a 
time-averaged  problem  was  basically  justified.  Time  dependent  effects  were 
essentially  of  two  types:  a  small  (0(cq)),  rapidly  oscillating  correction  to 
the  amplitude  and  an  0(e2)  correction  to  the  frequency.  This  analysis  breaks 


We  showed  how  to  treat  a 


down  as  q  or  6  becomes  comparable  to  e” 1 . 
representative  case  using  a  semiclassical  approximation.  It  was  found  that 
again  the  behavior  is  quasiperiodic  but  under  certain  "resonance"  conditions 
focusing  leading  to  large  localized  maxima  of  the  probability  amplitude  is 
possible. 

One  obvious  advantage  of  an  analysis  in  the  weak  coupling  regime  is  that  it 
provides  a  solution  that  can  easily  be  used  to  study  scattering  probability  as  a 
function  of  any  of  the  parameters.  For  instance,  we  provide  the  first  study  of 
the  sensitivity  of  scattering  to  the  standing  wave  intensity.  From  our  studies, 
we  find  the  inconclusive  results  of  the  experiments  not  to  be  surprising. 

Finally,  we  would  like  to  point  out  the  general  applicability  of  the 
perturbation  technique  employed  in  our  study  to  other  nonseparable 
time-dependent  situations.  In  fact,  the  qualitative  aspects  of  our  results  only 
depend  on  the  potential  being  a  periodic  function  of  space  and  time.  Extension 
to  more  general  time  periodic  potentials  is  straightforward.  Of  course,  results 
of  this  type  are  subject,  especially  in  their  long-time  applicability,  to  the 
limitations  of  multiple-scaling  perturbation  theory. 
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A^schematic  of  Kapitza-Dirac  scattering. 

Scattering  probability  P  as  a  function  of  q  when  t  =  2  for  scattering 
form  pi  =  hk  to  (a)  pf  =  3hk,  (b)  pf  =  hk,  (c)  *  -hk,  (d) 

pf  =  -3hk. 

Scattering  probability  P  as  a  function  of  q  when  t  *  2  for  scattering 
from  p£  =  -3hk  to  3hk. 

Scattering  probability  P  as  a  function  of  t  for  scattering  from 
p-  ■  -hk  to  pf  =*  hk  when  (a)  q  *  3,  (b)  q  *  5. 

Scattering  probability  P  as  a  function  of  t  for  scattering  from 
Pj^  *  -3hk  to  pf  =  3hk  when  (a)  q  »  3,  (b)  q  =  5. 

Scattering  probability  P  as  a  function  of  6  for  the  transition  from 
Pj  «  Bhk  to  =  (8  -  2)hk  when  t  =  2  and  (a)  q  *  3,  (b)  q  *  6. 

(a)  Phase  plane  for  Eq.  (4.12) 

(b)  Rays  (Eq.  4.11a)  in  x  -  t  plane 

The  surface  p  *  p(x,t)  (scales  as  shown). 

Structure  of  solution  near  a  fold  of  the  surface  p  =  p(x,t) . 

(a)  Outer  (classically  accessible)  region:  solution  can  be  given  as 
tH&'  superposition  of  two  waves. 

(b)  Boundary  layer:  solution  is  an  Airy  function  oscillating  with 
large  amplitude  to  the  left  of  X£  (caustic)  and  decaying  to  the  right, 

(c)  Outer  (classically  inaccessible)  region:  solution  decays 
exponentially. 

This  is  only  a  schematic  true  A  complex. 

Ray-caustic  coordinates:  r  is  arc  length  along  caustic  (solid  line), 
s  is  arclength  along  ray  (dashed  line).  At  C,  s  *  0  while  at  A, 
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